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processes the knowledge of the hadron distribution amplitudes, which are universal for different 
reactions, is essential. Only rather indirect information on these nonperturbative functions can be 
obtained from measurements. In this work we report on a lattice QCD computation of moments 
of nucleon distribution amplitudes using suitable three-quark operators. However, these operators 
have to be renormalized and the mixing is even more complicated than in the continuum. Using 
the symmetry group of the hypercubic lattice we therefore derive and implement irreducibly trans- 
forming three-quark operators, which allow us to control the mixing pattern and will finally lead 
to quantitative predictions in the MS scheme. We present preliminary results for leading-twist and 
next-to-leading twist nucleon distribution amplitudes based on the QCDSF/UKQCD simulations 
with 2 flavors of dynamical clover fermions. 
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1. Introduction 

QCD is the theory of the strong interaction, but nevertheless many results are obtained from 
perturbative calculations applicable to a large variety of scattering processes. The success of this 
approach is based on the factorization properties of the investigated reactions, which allow one to 
introduce distribution and fragmentation functions for quarks and gluons in the case of inclusive, 
and quark distribution amplitudes in the case of hard exclusive processes. These universal functions 
are common to different processes in all orders of the perturbative expansion. In ongoing and 
future experiments the investigation of hard exclusive and semi-exclusive processes will become 
the key tool for increasing our knowledge and understanding of the internal and spin structure of 
hadrons. Thus for facilitating fully quantitative predictions for these processes the knowledge of 
the nonperturbative quark distribution amplitudes is essential. They describe the hadron structure in 
terms of valence quark Fock states at small transverse separation and, unlike distribution functions 
in inclusive processes, cannot be accessed directly in experiment. Only some indirect insight can 
be obtained by measuring physical quantities like the magnetic form factor of the nucleon Gm{Q^)- 
At very large values of the electromagnetic form factors of the nucleon can be expressed as a 
convolution of the hard scattering kernel h{xi,yi,Q^) and the quark distribution amplitude in the 
nucleon (p{xi,Q^) [Q]: 

Gm{Q') = f\dx] f\dy](p*{yi,Q^)h{xi,yi,Q^)(p{xi,Q^) + 0{m^/Q^), (1.1) 
Jo Jo 

where [dx] = djcidjC2dx35(l — LLi equals the modulus of the squared momentum trans- 

fer in the hard scattering process. In this case only the leading twist nucleon distribution amplitude 
contributes. In an appropriate gauge, x, (yt) can be interpreted as the momentum fractions carried 
by the valence quarks before (after) the hard scattering. 

Apart from QCD sum rule determinations, an analytic approach to the distribution ampli- 
tude is feasible only for sufficiently large values of Q^, where the asymptotic form (p{Q^ ^ oo) = 
120x1X2X3 [§, ^ is obtained. However, given the logarithmic evolution in this knowledge is not 
really useful at reasonable energy scales, such that a nonperturbative lattice calculation seems to be 
the method of choice. 

At intermediate values of the momentum transfer (l < < lOGeV^) the electromagnetic 
form factors can be calculated from the nucleon distribution amplitudes using lightcone sum rules. 
In this case also higher twist terms of the nucleon distribution amplitudes ^ ^ will contribute. 
This kinematic region gained a lot of interest in recent years, because new data from JLAB |j8|, ^, [T^, 



11] for the well-known electromagnetic form factors of the nucleon contradict common textbook 



knowledge, for details see [12] and references therein. 

In this paper we want to present the theoretical framework needed to set up a calculation of the 
nucleon distribution amplitudes on the lattice and present some prehminary results for the leading 
and some next-to-leading twist distribution amplitudes. 

2. Theoretical background 

We would like to stress once more that the nonperturbative nucleon distribution amplitudes 
are needed anytime a process is calculated in which three quarks form a nucleon, e.g., after a hard 
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scattering process as displayed in Fig. |TJ In Minkowski space our starting point for the derivation 
of quark distribution amplitudes is the matrix element of a tri-local operator, 



(0| 



xd'y{z3)\Nip))e 



abc 



(2.1) 

where path ordering is implied for the exponentials, a,b,c are the color indices and \N{p)) de- 
notes the nucleon state with momentum p. We will consider these matrix elements for space-time 
separations of the quarks on the light cone with zt = aiz (z^ = 0) and = 1. 

Using the transformation properties of the fields in eq. ( ^ ) under Lorentz symmetry and 
parity it is possible to rewrite the leading twist contribution in terms of three invariant functions V, 
A and T [pi], 



4/iV { {p ■ YC)ap {Y5N)yVizi ■p) + {p- yy5C)aiiNyA{zi ■ p) 
+{i(y^vP''C)ap{f75N)yT{zi-p)}, 



(2.2) 



where C is the charge conjugation matrix, the nucleon spinor and fy the nucleon decay constant. 
Beyond the large limit there are also higher twist contributions, which we do not discuss here. 
In momentum space 



V{Xi) = V{xi,X2,X3) 



(2.3) 



the distribution amplitudes V{xi), A(x,) and r(x,) describe the quark distribution inside the nucleon 
as a function of the longitudinal momentum fractions x,. To be precise we want to note that the 
distribution amplitudes also depend on the factorization and renormalization scales. Here it is 
sufficient to set both scales to the same value jJ.. 
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Since so far it is not possible to access the quark distribution amplitudes directly on the lattice 
we consider their moments, which are defined as 

V''™= f\dx]x\4'x"^V{xi,X2,X3) (2.4) 

JO 



with equivalent definitions for the other distribution amplitudes. Using eq. ( [2.2| ) and ( |2.3| ) one can 
relate these moments of the quark distribution amplitudes to matrix elements of the following local 
operators 

r/''""(o) = £"^^'[^'i^^^ . .^^'M^(o)](c/)„/3r^^^ . .^^"'4(0)] 

X [/"^^' ...^^'■(75J''(0))t] (2.5) 

X [/"^^' ...^^"4(0)] (2.6) 
^P''^"{0) = e''^'\i'2^' . . . ^^'<(0)] ( C{-io^P)] [/'"^^' . . . ^^"'M^(0)] 

X [/"^^'...^^"(r^rsJ'-co)),] (2.7) 

by 

{0\rf'^mN{p)) = -fNV"""pPp'p'''p'N,{p), (2.8) 

{Ol^f'^^'mNip)) = -fyA''""pPp^p''p"N,{p), (2.9) 

{0\^f'^'mN{p)) = 2fyT"""pPp'p'"p"N,{p). (2.10) 

The multiindex Imn denotes the Lorentz structure given by the covariant derivatives = — 
igA^ in the local operator while p and T belong to a gamma matrix and a spinor respectively. The 
nucleon decay constant is normalized here by the choice V""" = 1 . 

For a lattice calculation the above equations and operators are converted from Minkowski to 
Euclidean space-time, then the techniques known from meson and hadron spectroscopy can be 
used to determine the matrix elements. However in the discretized space-time we are faced with 
additional obstacles being absent in continuum calculations. Though not mentioned explicitly, the 
above equations contain renormalization factors and the results obtained on the lattice should be 
converted to an appropriate continuum renormalization scheme such as MS. Due to the reduced 
symmetry of the discretized space-time we will have to treat additional mixings of operators with 
the same mass dimension which are forbidden in the continuum. Even worse, the discretized theory 
furthermore suffers from potential mixing with lower-dimensional operators. So it is important to 
notice that there is some freedom to choose especially well-suited operators exploiting the fact that 



different operators with distinct multiindices on the left hand side of eqs. ( |2.8[ )-( [2.10| ) are related to 
the same moments on the right hand side. Hence the first task is to find appropriate combinations 
of three-quark operators, for which the unwanted mixing is absent or at least strongly reduced. 

3. Irreducible representations of three-quark operators 



The matrix elements of the type {0\^\N) (p^-( O0| ) considered in the previous section refer 



to quark distribution amplitudes for the nucleon. However we want to stress that the results derived 
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in this section are more general and can be used in any kind of calculation involving local three- 
quark operators of the form 

^^''H^) = ^i2rM,..M„"«W"/3W^A'. •••^M»^rW' (3-1) 

where T^'^ is a tensor that represents the appropriate coefficients. As the actual position of the 
covariant derivatives does not influence the following discussion, we will assume that they act on 
the last quark unless stated otherwise. Isospin symmetrization will be discussed in detail in the 
following section. 

These three-quark operators are subject to renormalization and possible mixing under renor- 
malization. In order to get quantitative results from lattice simulations it is hence essential to 
perform a detailed study of their renormalization and mixing coefficients, preferably in a nonper- 
turbative approach. Let i^^^^ denote the lattice regularized, bare three-quark operator and its 
renormalized counterpart. They are related by a renormalization matrix Z: 

^(')-'-en =Z,-^-^W'''^^ (3.2) 

Operator mixing shows up in non-vanishing off-diagonal elements of Z. Typically several hundred 
independent operators may appear on the right hand side so that an elaborate approach is 



needed to gain control of the mixing issue. As in the case of quark-antiquark operators [ |15| ] the 
symmetry group H{A) of the hypercubic lattice provides appropriate tools to reduce the dimen- 
sion of the problem: One decomposes the operator space into subspaces transforming irreducibly 
with respect to the hypercubic group. Mixing under renormalization is then possible exclusively 
between equivalent irreducible representations, i.e., between operators that obey exactly the same 
transformation laws under the group action. Thus the Z-matrix becomes block diagonal and one 
has to care only about the lower-dimensional non-identical blocks, which typically mix an order of 
one to ten irreducible operators. 

We employ a similar group-theoretical approach for our three-quark operators. However, since 
half integer spin is assigned to our operators, we have to use the double cover of H{A), the so-called 
spinorial hypercubic group H{A) which was studied first by Dai and Song in 2001 This finite 
group contains 768 elements and can be defined by six generators, t, y, I\, h, h and U, and a set of 
generating equations: 



If 


= -1, 




= -ijii 


th 


= ht, 


th 


= ht, 




= ht, 


th 


= ht, 


rh 


= -/3, 


yh 


= -hr, 


rh 


= -hr, 


f 


= -1, 




= -1, 


{tyf 


= -1. 



Given these relations, it can be shown that, apart from the already known representations that are 
inherited from H{4) |[T8|], exactly five further inequivalent irreducible representations exist: zf, T^, 

Tp and T^. The superscript denotes their dimension and the subscript numbers inequivalent 
irreducible representations of the same dimension. 

While the representation describes how spinors and thus quark fields are transforming un- 
der group action, Lorentz vectors such as derivatives transform according to the representation 
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inherited from H{A). Thus one might in principle construct the transformation matrices for three- 
quark operators under action of any given group element G: 

^(;), transformed _ Q..ff(i)_ (3 3) 

However, this would again yield matrices of quite unhandy dimension (remember the amount of 
independent operators involved), so in a first step we reduced the set of possibly mixing operators 
and thereby the dimension of the transformation matrices by a detour via SO4. Writing the quark 
fields with dotted and undotted indices in the chiral Weyl representation |[l^] we can construct 
irreducible representations due to the local homomorphism SU (2) x SU (2) ~ S0{4). Each quark 
spinor naturally consists of two two-component spinors, and after contraction with Pauli matrices 
also the covariant derivatives can be written in an SU{2) representation: —>■ {^^o^Y^^. For 
any given operator, independent symmetrization of dotted and undotted indices projects onto an 
irreducibly transforming leading twist operator, e.g.: 

Uau''%d' Uau\&0)\d' U{aU^\&oY ^}d'\ (3.4) 

Fixing the indices in the above example one can read off twelve independent operators belonging 
to one 5C?4-irreducible representation. The chirality partners with exchanged dotted and undotted 
indices on the quark fields are treated in the same way. 

It is now straightforward to construct irreducible representations for the symmetry group O4 of 
the Euclidean continuum. As O4 = SOi[ U rSO/[ with r representing some reflection operation [jl^], 
joining the parity partners of the 5O4 irreducible multiplets into one generates an O4 irreducibly 
transforming multiplet. 

At that point the number of possibly mixing operators is sufficiently reduced to construct their 
768 transformation matrices Gij in eq. ( |3.3| ). Forming suitable linear combinations with the help of 
the characters of the spinorial hypercubic group H{4) we construct a projector P" that, applied to 
each O4 irreducible operator multiplet, projects out such operators that are in fact // (4) -irreducible: 

^" = ^ I r(Gr-G. (3.5) 

Here da denotes the dimension of the // (4) -irreducible representation to be projected on, ;t;"(G) 
the character of the group element G in that representation and |//(4)| = 768 is the group order. 
We present two typical examples for H{A) irreducible operators. First we display a subset of a 

12 

Tj— irreducible multiplet of operators with one derivative: 



2^2 ' ' 2V2 

^,u,uH^oy,,d^^-^^ 
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(mass)dimension 9/2 


dimension 11/2 


dimension 13/2 




derivatives 


1 derivative 


2 derivatives 


1 


5 multiplets 




3 multiplets 


4 






3 multiplets 




1 multiplet 


1 multiplet 


3 multiplets 




3 multiplets 


3 multiplets 


4 multiplets 






4 multiplets 


5 multiplets 



Table 1: Number of equivalent H{A) irreducible multiplets for three-quark operators with zero to two deriva- 
tives (not isospin symmetrized, derivatives acting on the last quark only). 



The curly brackets denote total symmetrization in the (un)dotted indices. As a second example we 

12 

show operators with two derivatives from one of the T— representations: 



Ordering the multiplets of irreducibly transforming three-quark operators according to their 
mass-dimension (number of derivatives) yields table This facilitates reading off the important 
results because mixing is possible only within one representation. As already mentioned, on the lat- 
tice operators of the same representation may additionally mix with those of lower mass dimension, 
e.g. involving an extra factor of l/a: 

^((),ren 2. .ffij)^^^'^^ _|_ 2' • — • ff^^^^ lower dim ^ g-^ 

However it is difficult and numerically challenging to factor out the divergent part. Therefore such 
operators should be avoided wherever possible or otherwise studied in detail so that the divergence 
can be cleanly isolated and finally subtracted from the renormalized operator. 

In our approach to the nucleon distribution amplitudes we exclusively use the representations 
Tp for three-quark operators without derivative, tj^ and for three-quark operators with one and 
two derivatives, respectively. Thus we safely circumvent any power divergences in terms of the 
lattice spacing. 



4. Irreducible representations of the nucleon distribution amplitude operators 

In the last section three-quark operators which are suitable for the calculation of matrix ele- 
ments on the lattice were derived. The next task is to establish a connection between these three- 
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quark operators and the quark distribution operators we are actually interested in. Before we give 
the corresponding relations, let us consider the additional symmetries among the moments of the 
quark distribution amplitudes introduced by the presence of two M-quarks in the nucleon: 



If we define 



which is a natural combination in our analysis, then the fact that the nucleon has isospin 1/2 implies 

r'™ = + (4.3) 

6 

With the analogous identities for V and A we can express the moments of V , A and T in terms of 
only one independent distribution amplitude 

ylmn ^1 (^^Imn ^ 20'"'" _ 0"''" _ ^"""'^ ^ (4.4) 
j^lmn ^1 (_2^lmn ^ 20™'" _ + ^"'A . (4.5) 



The combination' (p'™" = V'""' — A''"", often used in QCD sum rule calculations, can easily be 
expressed in terms of 0'™": 

Due to momentum conservation there exists also a connection between lower and higher moments, 

^Imn _ 0(/+l)mn _|_ ^l{m+l)n _|_ ^lm(n+\) ^ ^4^-7^ 

allowing us to test our calculation. 

Finally we relate the irreducible three-quark operators to the operators for the moments of the 
distribution amplitudes obtaining, e.g., 

= 1 (_^/3 ^ -^14 ^ .^23 ^ ^24 _ ^31 ^ -^32 ^ -^41 ^ ^42 (4.8) 

8 

where the lower index in the nucleon distribution amplitude operators denotes the spinor index and 

upper indices are Lorentz indices (see eq. ( |2.5[ )-(p7|)). The symmetry in the Lorentz indices on 

12 

the right-hand side reflects the leading twist projection of the operator. It is an element of the x— 

representation and hence does not mix with any operators of lower dimension. To construct an 

12 

isospin 1/2 operator we use an operator from another, equivalent x— representation, 

e^^^- =le"''' (-u\u\{^3dl)+iu\u\{&^dl)+u\u\{&xdl)-iu\^^^^ 

1 ^ ' (4.9) 

= - _ ^-^-M _ . ^^23 _ ^^24 ^ ^^31 _ • ^^32 _ -^41 _ ^42) ^ 



'Note that the normahsation conditions of (p"^" = 1 and (|)^"" = 3 are different 
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Combining both operators we obtain an operator of isospin 1/2 with its matrix element given by 

-4(0Ki°"+ - = A^i (pi - ip2) {P3 - ipA)fN (yoio - aOIO + 2r00^) . (4.10) 

In the same way we deduced the following set of operators used in our analysis of the nucleon 
distribution amplitudes moments on the lattice 



12 

• derivatives (Tp) 



ml^\N{p))=fy{piyi -P2n)N{p), 

m^i'mp)) = fN{P3n-pm)N{p), 

{o\^^^'^\N{p),) = Mpiyi +P2r2-P3r3-PAr4)N{p), 



12 

• 1 derivative (t— ) 



{0\i^l^\N{p)) =/w0"'" [{Yipi - r2P2)iY3P3 + Y4Pa) - 2pip2riYiMp), 

(OK|Af(p)) =/w0"'" [{7lPl + Y2P2){r3P3 - Y4PA)-2p3P4m4]N{p), 

ml^''\N{p)) =fN<l>''^{riPi - np2){nP3 - Y4P4)n{p), 



• 2 derivatives (t^) 

{01^2^^"^ \N{p)) = fN(l>^^^ {piP27l72 {P373 -P474) + P3P47374 {Pl7\ -P272)]N{p) 

with analogous expressions for other moments. The operators ffk denote isospin 1/2 combinations 
with k derivatives carrying an implicit spinor index. For the 0th moment we use only those operators 
which do not require non-zero spatial momenta. In the case of higher moments the advantage 
of using irreducible representations is paid by the requirement of non-zero spatial momenta. In 
the case of one derivative we used pi / and for two derivatives P2,P3 7^ 0. It is important to 
notice that for the 0th and 1st moments we can use several operators thus improving our statistics. 
Unfortunately for the 2nd moment we have to use operators from the representation to avoid 
mixing with lower-dimensional operators and hence we cannot increase our statistics in the same 
way. 

5. Unrenormalized results for the nucleon distribution amplitudes 

After these preparatory works we are finally able to obtain the desired moments of nucleon 
distribution amplitudes. We performed a lattice calculating of the above matrix elements, using 
parity and momentum projection for the appropriately chosen spatial momentum combinations. 
For the interpolating operator of the nucleon field we chose ^ = e"^'^ {u''CY5d^)u^^ leading to e.g. 

()ir4<M^(0,i;)>-/«»'"v/W "'''^'''''""+/^if-^"^ + ''^> exp[-£(S.] (5.1) 

where parity projection is implied. The normalization constant Zn{p) can be extracted from the 
usual nucleon 2-point correlation function 

{^{t,p)J^{0,p)) =Z^(p)^^^^^g^exp[-£(p>]. (5.2) 



9 



Moments ofnucleon distribution amplitudes 



T. Kaltenbrunner, N. Warkentin 



0.03 r 



> 

O 0.015 



■ p=5.29 ( 16" X 32) 

f^ = 0.00592 + 9e-05 GeV' 

♦ p=5.29 (24'x48) 



0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 

m " [Gev'l 




• p= 5.29 ( 16 X 32 ) 
-X| = 0.037+ 0.002 GeV^ 

* P= 5.29 ( 24' X 48 ) 
■ p=5.29( 16'x32) 

A, = 0.075 ± 0.004 Gev' 

♦ p= 5.29 { 24' X 48 ) 



J , , L 



0,3 0,6 0,9 1,2 1,5 1,8 2,1 2,4 

m " [GeV^] 



Figure 2: Nucleon decay constant fy (upper plot) and the normalization constants —Xi and X2 of the 
next-to-leading twist distribution ampUdtudes (lower plot). All displayed results are unrenormalized. 



Since our operators for the moments of the nucleon distribution amplitudes are local we apply 
smearing only to the nucleon interpolating operator at the source. The two-point nucleon correlator 
is smeared both at source and sink. Due to this procedure the location of the effective mass plateaus 
is different for the two-point nucleon correlator and the two-point correlators used to determine the 
moments. Hence we do not calculate the ratios of these operators, but perform fully correlated 
fits to the correlators themselves choosing the fit ranges to match the effective mass plateaus. We 
use the pion masses determined by the QCDSF collaboration and set the scale with the Sommer 
parameter tq = 0.467fm. 

For the sake of flexibility we perform a two stage analysis. In the first step we calculate two- 
point functions of general 3-quark operators schematically given by 

(0|£''^''[^^' ...^^"J]"^^r(;7)|0), (5.3) 

with J + m + n <2. In a second step, these general operators can be used to calculate different 
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■ P= 5.29 ( 24' X 48 ) 

(])"'=0.870+ 0.006 

♦ P = 5.29 (24"' X 48) 



<> 



, , , , , , , , , , , I , , \ , , I , , I , , L 

0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 

ra^ [Gev'] 



0.5 
0.4 
0.3 


1 1 1 1 1 , 1 1 1 1 1 , 1 1 1 1 1 , 1 1 , 1 1 

■ |i=5.29( 16"'x32) 
(|>™ = 0.27 ± 0.02 


♦ p= 5.29 ( 24' X 48 ) 












„ -\ 


...... 




0.2 
0.1 










. . 1 . . 1 . . 1 . . 1 . . 1 . . 1 . . 1 . . 



m^' [Gev'] 

Figure 3: 0""' (upper plot) and 0"^" (lower plot) moments of the leading twist nucleon dsitribution ampli- 
tude. All displayed results are unrenormalized. 



1.1 
1 

0.9 
0.8 
0.7 



matrix elements, in our case the irreducible combinations relevant for the moments of nucleon 
distribution amplitudes. 

Although not discussed here, the calculation of moments of higher twist distribution ampli- 
tudes follows the same procedure as for leading twist. Thus using the general three-quark operators 
from step one, we can extract, e.g., the normalization constants X\ and A2 of the next-to-leading 
twist distribution amplitudes, which are also related to matrix elements of local three-quark opera- 
tors [|]. 

Some of our results along with simple linear chiral extrapolations are presented in Figs. ^ and|3[ 
In Fig. |2| (upper plot) we show the normalization constant of the nucleon distribution amplitude as 
function of the pion mass. The analysis was done on two different volumes, 16^ x 32 and 24^ x 48. 
We observe small finite size effects on the 16^ x 32 lattice for the smallest pion mass. In the lower 
plot of Fig. ^ we give —X\ and A2, the normalization constants of the next-to-leading twist distri- 
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bution amplitudes. Compared to the leading twist case, we observe more pronounced finite size 
effects for these quanitites. 

In Fig. ^ we give representative results for two of the higher moments, <^^^ (upper plot) and 
^020 (lower plot), as a function of the pion mass. The data for the 16^ x 32 and 24^ x 48 lattices 
do not seem to be completely consistent. Improving statistics, we are currently examining this 
behavior, which is observed only for some of the higher moments. 

The results obtained on the 24-' x 48 lattices are extrapolated linearly to the chiral limit. The 
resulting values for the normalization constants fy, X\ and A2 are close to the QCD sum rules 



calculations [19, EOl], but after renormalization we expect them to be lowered by approximately 



twenty percent. 



6. Renormalization 

All results presented so far refer to unrenormalized lattice calculations. The important step 
of renormalization for the three-quark operators used is in progress and will be discussed in a 
forthcoming paper, while here we will only sketch the keystones. For the renormalization we 
adopt a nonperturbative procedure, analogous to the well-established RI —MOM scheme Martinelli 
et al. introduced for quark-antiquark operators [^. We contract our isospin symmetrized and 
color antisymmetrized three-quark operators with three quark momentum sources on the lattice 
and calculate a correlation function of the following kind: 

G^'\p,q,ry^'^^ = ldudvdwdxe''^'"+P''+''''^e-'^P+''+'>'{u{u^^^ (6.1) 

with a, P and / denoting spinor indices and p, q and r being the incoming quark momenta. The 
gauge is fixed to Landau gauge. After amputating the external legs, we are left with a three-quark 
vertex F that contains all radiative corrections: 

G« {p,q, r)%^ = F« {p,q, r)„,pyS{-r)a,ccS{-p)p>pS{-q)y,e'"'' . (6.2) 

Imposing a suitable renormalization condition at the scale jx^ = (p^ +q^ + r^)/3 we obtain the 
renormahzed vertex 

F«'-"(;.,^,r;M)=zJ(At)-Z,'/'(Ai)-FW(/.,^,r). (6.3) 

The RI — MOM renormalization matrix is then converted to the MS scheme using one loop 
continuum perturbation theory. Finally we can perform a renormalization group extrapolation to 
any renormalization scale jJ. we like. 

As an outlook and in order to demonstrate how well this procedure works, we want to present 
a consistency check for the zeroth and first moments of our nucleon distribution amplitudes. Due 
to momentum conservation the sum of the three first moments is equal to the zeroth moment (cf. 
(0)): 

0100 + ^010^^001^^000 ^3^ (6.4) 

One can add up the bare first lattice moments and will receive a value of 2.60, which is incompatible 
with the above sum rule (fig. 0, dashed line). Incorporating the three-quark operator renormaliza- 
tion and especially the mixing matrix elements (we have to deal with a 3 x 3 matrix in this case) 
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Figure 4: Consistency check for the zeroth and first moments. Dashed line: sum of the bare first moments. 
Full line: renormalized sum of the first moments. The moments renormalized at different scales fit the 
theoretical constraint very well. 



we obtain very good matching. We would like to emphasize that the well reproduced mean value 
of 3.00 results from a renormalization at different scales and that the error bars containing uncer- 
tainties in the ^s as well as in the Zs are a rather conservative estimate at this stage. This check 
makes us confident that our preliminary results tend into the right direction and deserve further 
investigations. 



7. Conclusions and perspectives 

In this paper we have discussed how moments of nucleon distribution amplitudes can be cal- 
culated on the lattice. After rewriting the lightcone expressions in terms of local matrix elements, 
we focused on operators that are well-behaved even in the context of reduced lattice symmetries. 
The representations of the spinorial hypercubic group provide the tool to derive irreducible sets 
of three-quark operators that are suited to control mixing and renormalization. After appropriate 
isospin symmetrization we expressed the operators for the moments of the nucleon distribution 
amplitudes in terms of these irreducible representations and showed preliminary unrenormalized 
results for various moments of the distribution amplitudes. The renormalized results for the first 
moments are in good agreement with a sum rule. A comprehensive and detailed study of the three- 
quark operator renormalization is in progress right now. 
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